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\Q | Abstract 

The spectrum of electromagnetic fields satisfying perfectly conducting boundary conditions in a segment 
i | , of a straight beam pipe with a circular cross-section is discussed as a function of various source models. 

■ These include charged bunches that move along the axis of the pipe with constant speed for which an exact 
£Lij solution to the initial-boundary value problem for Maxwell's equations in the beam pipe is derived. In 

■ the ultra-relativistic limit all longitudinal components of the fields tend to zero and the spectral content 
' of the transverse fields and average total electromagnetic energy crossing any section of the beam pipe are 

directly related to the properties of the ultra-relativistic source. It is shown that for axially symmetric 
ryj ■ ultra-relativistic bunches interference effects occur that show a striking resemblance to those that occur 

O ' due to CSR in cyclic machines despite the fact that in this limit the source is no longer accelerating. 

' ^2 , The results offer an analytic description showing how such enhanced spectral behaviour depends on the 

geometry of the source and the details of the stochastic distribution of structure within the source. The 
' field energy spectra associated with a source containing Af identically charged ultra-relativistic pulses, each 

Oh with individual longitudinal gaussian profiles distributed according to a uniform probability distribution 

with compact support, is compared with that generated by charged bunches containing a distribution with 
I , 2n + 1 peaks in a region with compact support (modelling micro-bunches). These results offer a viable 

■ experimental means for inferring properties of the longitudinal charge distribution of bunches with micro- 
' structure in ultra-relativistic motion in straight segments of a beam pipe from observation of the associated 

ff) , electromagnetic energy spectra. They are also of relevance to design criteria where the coherence effects of 

OO ■ such fields play a significant role. 

cn ; 
od ■ 

o : 

1 Introduction 



PACS 87.56.bd, 02.40.-k, 41.20.-q, 29.27.-a, 41.60.-m, 41.75.-i 



k>( | Modeling the behaviour of charged particles in a modern accelerator is a critical component in its design. Due 
; i ■ to the inherent non-linear nature of the dynamics of charged bunches in external electromagnetic fields such 
! modeling generally necessitates numerical computation. With current hardware such computations often require 
further stringent approximations on the equations of motion in order to extract viable information. To date 
the most sophisticated Maxwell- Vlasov solvers are unable to effectively model electromagnetic interference in 
3 spatial dimensions. Approximations often neglect the effects of confining boundaries and wakes, radiation 
reaction forces, detailed stochastic properties of the beam and possible quantum effects. However, in certain 
limits analytic information can be deduced from the fundamental equations of motion for the coupled particle 
field system. If the motion is prescribed the problem reduces to solving Maxwell's equations for convected 
sources. Nevertheless finding solutions satisfying boundary condition appropriate to an accelerator is in general 
a non-trivial exercise. 

Since Fourier analysis is a linear operation the spectral content of the electromagnetic energy in the fields 
produced by charged particles in prescribed motion in free space has been exhaustively investigated over many 
decades, particularly for collections of particles in uniform circular motion [UI2]. The effects of the superposition 
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of retarded free-space solutions to Maxwell's equations with distributed currents on the spectral content of their 
radiation fields gave rise to the notion of coherent synchrotron radiation (CSR) in which enhanced radiation 
in some frequency domain can occur, depending quadratically on the total number of point particles in an 
accelerating bunch [3] . It is generally assumed that the effects of confining boundaries in a real cyclic machine 
do not significantly alter the criteria for CSR from those that arise for sources in free space. However analytic 
efforts to determine such criteria taking account of boundary conditions and stochastic effects inevitably demand 
further approximations [H[S]. Since the role of CSR is a fundamental ingredient in the design of new light sources 
it is of interest to explore new analytic approximations schemes that can complement the numerical simulations 
of particle-field interactions. 

In this note we approach the modeling process by ensuring at the outset that all fields satisfy perfectly 
conducting boundary conditions in a segment of a straight beam pipe with a circular cross-section. The sources 
will be assumed to move along the axis of the pipe with constant speed and an exact solution to the initial- 
boundary value problem for Maxwell's equations in the beam pipe is derived. In the ultra-relativistic limit all 
longitudinal components of the fields tend to zero and the spectral content of the transverse fields and average 
total electromagnetic energy crossing any section of the beam pipe are directly related to the properties of the 
ultra-relativistic source. Thus a measurement of the former offers a direct method of estimating properties of 
the latter see e.g. [Hj. However the detailed structure of the spectral content depends on the source structure. In 
the following we demonstrate that for axially symmetric ultra-relativistic bunches interference effects occur that 
show a striking resemblance to those that occur due to CSR in cyclic machines despite the fact that in this limit 
the source is no longer accelerating. Furthermore it is possible to study analytically how such enhanced spectral 
behaviour depends on the geometry of the source and the details of the stochastic distribution of structure 
within the source. 

Section 2 establishes the formalism for solving Maxwell's equations for fields inside a perfectly conducting 
cylindrical beam pipe in terms of complex Dirichlct and Neumann eigen-modes of the 2-dimensional scalar 
Laplacian operator. An exact solution to the initial-boundary value problem is exhibited and this is used to 
explore fields associated with ultra-relativistic bunches. Section 3 develops the model in terms of an axially 
symmetric source, offers a compact formula for the spectral distribution of electromagnetic energy that crosses 
any section of the beam pipe and, in section 4, discusses the effects on the spectra produced by different types 
of stochastic distribution. 

2 Electromagnetic Field Solutions 

The electromagnetic fields e and h in a beam pipe in the presence of sources with charge density p and electric 
current density J satisfy the Maxwell system: 

V x e + [Lq d t h = 
V x h - fi Y 2 d t e - J = 
V • h = 
HoY 2 V • e - p = 

where the admittance Y = l/(poc) with c being the speed of light in vacuo. All vectors will be referred to 
a global ortho-normal frame {e r ,e z ,eg} defining a cylindrical coordinate system {r, 9, z} with the axis of the 
beam pipe along the z— axis. 

A transverse circular section of this pipe is denoted T> with boundary dT>. In the following we exploit 
the properties of a complex Dirichlet mode set {$at}. This is a collection of complex eigen-functions of the 2- 
dimcnsional (transverse) Laplacian operator V 2 on T> that vanishes on dT>. This boundary condition determines 
the associated (positive non-zero real) eigenvalues 0^. The label N here consists of an ordered pair of real 
numbers. Thus 

V2$ w _ /3 2 N $ N = 0, (1) 
with $iv|ax> = 0. With an overbar indicating complex conjugation these modes are normalised to satisfy 

/ Wi$Nrdrd0=Nl6 NM , (2) 



2 



An explicit form for is for neZ 

*N(r,0) = Jn (x q(n) -^ e m0 , (3) 

where J n (x) is the n-th order Bessel function and the numbers {xg( n )} arc dehncd by J n (xq(n)) — and 
N := {n, q(n)}. The eigenvalues are given by {fi^ = x q ^/a} and Af^ = ira 2 J 2 +1 (x q ^). 

In a similar manner a Neumann mode set {^n} is a collection of eigen- functions of the Laplacian operator V 2 
on T> such that ^f- vanishes on dV. This alternative boundary condition determines the associated (positive 
non-zero real) eigenvalues a 2 N where again the label N consists of an ordered pair of real numbers: 

V 2 *iv - = 0, (4) 

with ^f-\ov = 0. These modes are normalized to satisfy 

W^^ N rdr d6 = M 2 N S NM - 
An explicit form for ^> m is for m G Z 

<Mr,0) = J m (^ (m) ^)e™ e , (5) 

where the numbers {^p( m )} are defined by Jm( x 'p( m )) = an( l A/ := {m,p(m)}. The eigenvalues are given by 

{«« = 2>( m) / a } and >*M = TO2j m+iK( m ))- 

The electromagnetic fields in the interior of the cylindrical beam pipe of radius a satisfying perfectly con- 
ducting boundary conditions at r = a can now be expanded [SJ [S] as: 

e(t, z, r, 9) = Y, VnV, z) V$ N (r, 9) + £ y£(t, z) e z x W M (r, 0) + ^ 7*(*. *) $ M(r, 0) e, (6) 

JV M JV 



h(t, z, r,9) = Y I*(t, z) e z x V$iv(r, 0) + £ i£(t, z) V* M (r, 6) + £ 7 £(i, z) * M (r, 0) e 2 (7) 

TV M M 

The fields are assumed to be generated by an external RF source that accelerates charged bunches to near 
the speed of light. In a straight beam pipe such a source can be modeled by an arbitrary smooth convective 
charge density p and a current with ortho- normal components J r = Jg = 0, J z (z — vt,r,6) = v p(z — vt,r,9) 
with v close to the speed of light. The equations for 7^ and 7^ that follow from Maxwell's equations and (©, 
([7|) for these sources are: 

In - c 2 lT + c 2 a 2 Nl % = (8) 
2 

■• E 2 Ell , 2o2 E c MO / 2 2\~T~ ln\ 

7n- c 7n + c PnIn = — 772-(c - v )Pn (9) 

■/Vjy 

where 

p w := / p$ N rdrd9. (10) 
In terms of 7^ and 7^ and the projected convective sources 

(11) 
(12) 

(13) 
(14) 

Finding the fields e and h is now reduced to solving an initial- value problem for the decoupled fields 7^, 7^. 
For some real constant a > and source g(t, z) each is a solution to the generic (hyperbolic) partial differential 
equation 

j-c 2 .f" + c 2 a 2 f = g. (15) 



V N — 


1 1 

P 2 N { lN Af 2 N p 




v H - 


2 Tjv, 




1 N — 






T H 






1 N — 
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The general causal solution f(t, z) with prescribed values of /(0, z) and /(0, z), is (see e.g. [TU] ) 

/(t, 2) = W„ [/*""](*, 2) +2 a [ ff ](i, z), (16) 



where 



W CT [/'""*] (t, z):=Uf(0,z- ct) + /(0, z + ct) 



2 



1 

2c" 



^/(COJoK^-^-C) 2 ) 

2 — ci 



and 



1 



r*" J(/ (,,( l *^H!), (17) 

2 A-a v/?! 5 - (= - O 2 

t rz + c(t-t') 



2c 7 



l a [g}(t,z):=- dt' d(g(t'X)J (ay/cZ(t-t>)z~(z~O 2 l (18) 



z-c(t-t') 



The functions f(0,z),f(0,z) constitute the initial t = Cauchy data in this solution and determine the 7i a 
contribution above. 

For a bunch with total charge Q moving with speed v we assume here that p can be written 

p(z - vt, r,6) = Q p x (r, fl)pH (z - vt), 
where v is given (v < c) and /0 _L (r, 9), p"(z — ft) are arbitrary smooth functions subject to 

/oo 
dz p" (z — vt) = 1, 
-00 

With these sources the causal solutions to © and © for 7^ and 7^ are given by: 

7jV (M) = H aN [l§ mU ](t,z), 

~/§(t,z) = HpAlN nit }(t,z)-^(c 2 -v 2 )X 0N {^](t,z), 
where p N is given by (TIT)]) and 

p' N {t,z) = QpW'{z-vt) ( p^(r,6)$ N rdrde. (19) 



In the ultra-relativistic limit, v — ► c, the second term in 7j^(i, z) tends to zero. 

3 Spectral Energy Distributions 

When the transverse distribution depends only on r, expressions for the electromagnetic fields associated with 
the sources simplify. The source under consideration is axially symmetric if 

p ± (r,6)=K(r), (20) 

where lZ{r) is a smooth function satisfying 

r 1 

/ drrll{r) = —. (21) 
Jo 2?r 

Then for axially symmetric bunches 

p' N (t,z) = 2nQtot6n,oP l] '(z-vt) J dr rTZ{r)J (s 9 (o)~) ■ (22) 

From these formulae the instantaneous real electromagnetic power crossing any section of the beam pipe is 
obtained by integrating the component S = (e x h) • e z of the Poynting vector field over the cross-section T> at 
an arbitrary point with coordinate z 

V(t,z) := [ S(t,z,r,6)rdrd8 
Jv 
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With the aid of the orthogonality properties of the Dirichlet and Neumann modes this becomes 



V = k{ E^ F ^ -H< M « V m^\ (23) 

^ N M ' 



where 9? takes the real part of its argument. From (fTTj) . (fT2"|) . (|13p . and (fTl|) this power flux can be 

explicitly expressed in terms of the source projections pn, 7m an d In 



V = x(j2 

^ JV 



, El 



J^rflrltoY* + n Pn P n Pn 



"j&yS'tf) -moE#7M'|- (24) 

/ AT ) 



/ M "A/ 

In a cyclic machine with time-periodic fields of period T (such as a synchrotron) the fields associated with the 
sources can be expanded in a Fourier series: 

1 f T 

e n (z,r,0) = — / e(f, z, r, 6) exp(iQt) dt 
1 Jo 

1 f T 

h n (z, r, 9) = — / h(t, z, r, 9) exp(i£tt) dt 
T Jo 

where fi = 27r/T, and the mean (time-averaged) power crossing D at z is 

oo 

< v > o) = J2 < Vn > ( z ) 

n— — oo 

where 

<V n >(z)= I (e n (z,r,e) Xh n (z,r,9)) ■ e z )rdrdO 



\ > (z) = / (e n (z,r,9) x h n (z,r,9)) ■ e z )r 
Jv 



However in a straight beam pipe the fields associated with the sources are not periodic in time. The sources 
are however localized in space so the total electromagnetic energy hi crossing any section T> at z is well-defined: 

U(z)= V(t,z)dt= (—{oj,z))dcj 



doj 



where 

dll_ 

du) 

in terms of the Fourier transforms: 



(cj,z) = J (e((j,z,r,9)xh{u,z,r,9)-e z ^ rdrd9 (25) 



1 f°° 

e(w, z, r, 9) = —= I e(i, z, r, 9) exp 

V Zlt J-oo 



(iurt) dt 



h(cj, z, r, 9) = / h(t, z, r, 9) exp(iu;t) dt 



/27T 

Clearly ^ff (w, z) du; is the total electromagnetic energy crossing V at the station z of the beam pipe in the 
wave band w\ to u>2 and is experimentally accessible with suitable detector diagnostics. Measurement of the 
spectral content of such energy associated with ultra-relativistic charged bunches can offer valuable information 
about their longitudinal charge distribution. 

In the ultra-relativistic limit the fields e and h have components that lie solely in the transverse sections 
of the pipe and are concentrated in space where the sources are concentrated. For an axially symmetric ultra- 
relativistic bunch with charge Q one calculates from (f2"Sj) 

where the constant 

c 2V \f v K(r)<t> M (r,9)rdrdO\2 
6cT 4/ \^hi! a a J 2 {pMr)rdrW 
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In general one expects the source to require a stochastic description since its origin is fundamentally stochas- 
tic. This requires the introduction of stochastic variables and their associated probability distribution in order 
to calculate the expectation values of However it is of interest first to note the deterministic structure of 

that arises from a well-defined finite train of M charged ultra-relativistic pulses where each pulse has the 
same longitudinal profile / in p" . Denote the longitudinal charge distribution p" of such a train by F with 

M 

F(z-ct) = FoJ2f(z-c(t- T j)) 

3=1 

for some real constant Fq. 

If f(u>) is the Fourier transform of f(t) then the Fourier transform of Y^jLi /(* + T j) with respect to t is 
/MEjli exp(-iujTj) and 

N / N N \ 

|/H^cxpH^)| 2 -|/H| 2 hV + 22E«»w(7 1fe -r J (26) 

3=1 V 3=lfe=l / 



Hence for a train of such equidistant pulses with spatial separations cTq > 

^{u) = KF*\f\\u) C<u) 

where C(ui) = 1 \l°^uT\$ an( ^ lui ^ s independent of z. 

The function C is bounded (0 < C(u>) < TV 2 ) with maxima at u> = o>j = but the detailed behaviour of 



^ff(w) depends on the single pulse structure defined by /. Thus if 




describes the structure of the j — th pulse in the train with Tj = j Tq one finds 

dU — 2 



7T<7 



where 

1 - cos(wTqAA) 



W{uj) = exp - 




1 — cos(wTo) 



The modulation of C by |/| 2 means that the maxima of C |/| 2 are shifted from ui = ojj. For typical charged 
bunches described by the above Gaussian form for / the shift is unlikely to be experimentally detectable. 
However the points where u> = ujj determine the first order contact points of the curve ^ff(w) with the envelope 
curve 

£+H^A-F 2 |/| 2 HAA 2 
Lira -rn = 1 



dU 



Lim —wj- - - = 1 

The general features of the spectrum of ^ for the choice of Gaussian fj are sketched in figure 1. It should 
be stressed that the electromagnetic fields generated by the non-stochastic source under consideration arc fully 
"coherent" . The maxima in ^ are produced by constructive interference of the fields associated with the 
regular structure in the train that is maintained during its ultra-relativistic motion. The spacing of adjacent 
maxima of ^(w) produced by bunches in an accelerator differ imperceptibly from the spacing jr of adjacent 
points where ^ is tangent to £ + . Thus the longitudinal spatial separation cTq of the maxima in this idealized 
bunch containing a structure with Af equidistant peaks is immediately visible in the electromagnetic energy 
spectrum. 
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Electromagnetic Energy Spectrum 
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Figure 1: The rescaled electromagnetic energy spectrum with x = TqMld: 
W(x) = W( T x tf ) associated with the longitudinal source profile 



£iZa ^ \ _ p or purposes of visualization a z 



hcT , Af = 5, T = 1, 



21] 



The upper dotted curve denotes the envelope 



4 Stochastic Effects 

To see the effects of randomization on these interference maxima a simple stochastic model of a charged bunch 
containing Af identifiable random variables will be adopted. Instead of fixing the separation of the Af peaks 
in p" they will be distributed according to some probability measure . We choose as random variables 
T\ , T 2 , . . . T/v and assume that 



P N {T 1 ,T 2l ...T N ) = I^ =1 P{T J ) 

where P(T) is a probability distribution for a single random variable T. The expectation value of any function 
O of Ti,T 2 , . . .Ttf will be denoted E P u{0) where 

E P u{0) = f 0(T 1 ,T 2 ,...T^)P Af (T 1 ,T 2 ,...T^)dT 1 dT 2 ...dT Af 

In particular it follows from (|26[) that the expectation value of the spectral energy distribution is 

dU 



EpM ( ^(w) ) = A-F 2 |/| 2 H E p m{L(u)) 



E p h{C{lu))=N +{N 2 - TV) I / P(T)e^ T (ff| 



where 



The magnitude of (jV 2 — A/") times the modulus squared of the Fourier transform of P(T) in this expression 
relative to Af determines the nature of the expectation value of the spectral energy as a function of to. This 
expectation value is now bounded above and below by two distinct envelopes that vary with to. Within these 
envelopes one may classify local maxima as major and minor (see figure 2). The width of the first local dominant 
maxima in the u) spectrum of the expectation value is directly related to the overall scale of the spatial size 
of the bunch source as determined by the probability distribution P(T) while the separation between adjacent 
local major maxima in ^ is determined by the structure of P(T). These features are illustrated below where 
one notes that the bounding envelopes have single maxima at co = in the ratio Af 2 : AT. For large Af this is 
then close to the ratio of the first few ratios of ( local maxima : local minima) of §^ . Following tradition it 
is natural to refer to fields that contribute to the first few major local maxima of as exhibiting "stochastic 
coherence" . The relation between the non-zero frequency at which the first local minimum of ^ occurs (or the 
frequency beyond which ^ lies close to the lower bounding envelop) and the spatial distribution of charge in 
p" is a stochastic one depending on the structure of P(T). These general features are illustrated as follows. 
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• If Af identical pulses in p>>, each with the above longitudinal Gaussian profile /, are, for some constant F, 
independently distributed according to P N with: 

1 



P(T) = - for - - < T < - 
v i r 2 - - 2 



and zero elsewhere then 



|'* M )=^|/|V)U +( ^-A0(^ r/2 j 



iMV2)y 



The bounding envelopes (seefigure2) are the curves £ (w) = ifi* 1 , 2 A" |/| 2 (w) and£ + (u) = KF^Af 2 |/| 2 (w): 

£-(u) < EpM f^jjHj < 

It is natural to designate the Fourier components of fields as "stochastically incoherent" if they contribute 
to the expectation E p m (fff (w)) in the vicinity of the lower bounding envelope £~(ui). This expectation 
value first touches £~(uj) at approximately u = and V determines the average spatial length of p". 

Electromagnetic Energy Spectrum 
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Figure 2: The rescaled electromagnetic energy spectrum with x = Tlu: 

W{x) = E p h (jr)) /(FfjK) associated with the longitudinal probability distribution 

P(T) = 1/r for — T/2 <T< T/2. For purposes of visualization a z = ±cT, F = 1, c = 1, Af = 10. The upper 



dotted curve denotes the envelope fLiS-l and the lower one the envelope jwtj? 



• If, for some integer n and constant k, Af identical pulses in p" , each with the above longitudinal Gaussian 
profile /, are independently distributed according to P 1 ^ with: 



P(T) 



1 



(2n + 1)kJtt z 



I] exp 



T-jr 



for — oo < T < oo 



then 



= A"F 2 |/| a ( w ) ( + (AA 2 - AO 



(2n + l) 26XP I ^ 



1 - cos(wr(2n + 1)) 
1 — cos(W) 



Such a distribution P(T) offers a means to model the electromagnetic energy associated with a bunch 
containing Af constituents with 2n+l micro-bunches [11] distributed stochastically among them according 
to equally spaced Gaussian distributions with parameters k and r. 

In this case the bounding envelopes (see figure 3) are the curves 



£-{Lo)=KF 2 N\f\ 2 {u J ), 
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£+M = KFi |/| 2 (u,) [n-{N 2 -M)eM-^f)) 

and 

£~{w) < E p m f ^j(<*>)J < £ + M 

For w 7^ the envelope here lies below the envelope £ + (uj) for the previous distribution in figure 

2. Hence the overall effect of the micro-structure introduced in this P(T) is to suppress the magnitude of 
the j > local major maxima near u = in the expectation values of ^j. 



Electromagnetic Energy Spectrum 




x T 

Figure 3: The rescaled electromagnetic energy spectrum with x = kuj: W(x) = E p m ($rf- (-)) /(FgK) associ- 
ated with the probability distribution 
P(T) = ((2n + 1)«0F) E;=-„ exp (— (T - jr) 2 ^ 2 ) . 

For purposes of visualization cr z = | ck, r = 30K(n — , Af = 10, n = 2, k = 1, c = 1. The upper dotted 
curve denotes the envelope f F i^\ and the lower one the envelope fp&^X ■ 



5 Conclusions 

A general discussion of the electromagnetic energy spectra generated by ultra-relativistic distributions of charge 
in a perfectly conducting straight cylindrical beam pipe has been given. A detailed model has been constructed 
for particular distributions where the electromagnetic fields have been calculated ab-initio by solving analytically 
the initial-boundary value problem for Maxwell's equations in the beam pipe. The effects of the pipe geometry 
and the source proper charge distribution on the field structure have been explored for different stochastic models 
and the resemblance of this structure to CSR behaviour in a cyclic machine pointed out. Further aspects of 
such dynamically enhanced field effects can be found in |12j . Our results offer a viable experimental means for 
inferring properties of the longitudinal charge distribution of bunches with micro-structure in ultra-relativistic 
motion in straight segments of a beam pipe from observation of the associated electromagnetic energy spectra. 
They are also of relevance to design criteria where the coherence effects of such fields play a significant role. 
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